Intrinsic decoherence in atomic diffraction by standing light wave 
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We study the effects of time uncertainty in the interaction of atoms with a standing light wave. 
We discuss its physical origin and the possibility to observe intrinsic decoherence effects by measuring 
the atomic momentum distribution. 
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Typically the decoherence effects are ascribed to the action of an environment, i.e. not controlled degrees of freedom. 
This is the so called environment-induced decoherence Q. However, there could also be the possibility of an intrinsic 
O ji nature of decoherence, i.e. coming from the system itself and not from external degrees of freedom, as pointed out 
recently |Q. The idea underlying this approach is the fact that the interaction time, i.e. the time interval in which 
the effective Hamiltonian evolution takes place, is a random variable ||. This randomness can have different origins 
depending on the studied system. 
I , In atom optics Q], for instance, the interaction time is determined by the transit time of the atoms through a 
' cavity ||. Nevertheless, this interaction time is random, due to eventual classical fluctuations of atomic velocities, 
, but beyond that, due to the finite extend of the interacting atomic wave packet ||. In fact, the instant at which a 
wave packet passes a given point of the axis is not determined precisely, but carries an uncertainty proportional to 
its width Hence, there should be purely quantum mechanical intrinsic decoherence effects. Here, we study the 
appearance of such effects in atomic diffraction by standing light wave. 

On the other hand, decoherence effects relevant to quantum optics domain are now beginning to be tested exper- 
imentally 0| . As far as quantum objects bearing mass are concerned, decoherence has been only investigated for 
the motional states of the ions Then, it would be also interesting to experimentally verify such a new type of 
1 decoherence. To this end, we show a possible modification of the measurable atomic momentum distribution. 

CS ' II. TIME FORMALISM 

3 ' 

In this section we review the formalism describing non-dissipative decoherence derived in Refs. j|. It is based 
on the idea that time is a random variable or, alternatively, that the system Hamiltonian (therefore its eigenvalues) 
fluctuates. This leads to random phases in the energy eigenstates representation. Then, the resulting evolution of 
the system must be averaged on a suitable probability distribution, and this leads to the decay of the off-diagonal 
?3 , elements of the density operator. 

Let us consider an initial state p(0) and consider the case of a random evolution time. Then, the evolved state will 
be averaged over a probability distribution P(t,t'), i.e. 

poo 

p{t) = / dt'p(t,t')p(t'), (I) 



JO 

where p(t') = exp{— iLt'}p(0) is the usual solution of the Liouville-Von Neumann equation with L . . . — [H, . . ]/Ti. 
One can write as well 

p(t) = V L (t)p(0), (2) 

where 

POO 

V L (t) = / dt' P(t,t')e~ iLt ' . (3) 
Jo 

In Ref. ||, the function P(t,t') has been determined to satisfy the following conditions: i) p(t) must be a density 
operator, i.e. it must be self-adjoint, positive-definite, and with unit-trace. This leads to the condition that P(t,t') 
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must be non-negative and normalized, i.e. a probability density in t', so that Eq.([j]) is a completely positive map; ii) 
Vz,(f) satisfies the semigroup property Vr,(ti + t<z) — VL{ti)Vc(t2), with t\,t<i > 0. These requirements are satisfied by 

^=(1+^' (4) 



WHr^TTTZT r - (5) 



and 



TT(t/T) \t 

where the parameter r naturally appears as a scaling time. The expression (^|) is the so-called T-distribution function, 
well known in line theory. Its interpretation is particularly simple when t/r — r with r integer; in that case P(r,t') 
gives the probability density that the waiting time for r independent events is t' and r is the mean time interval 
between two events. Generally, the meaning of the parameter r can be understood by considering the mean of the 
evolution time (t 1 ) = t, and its variance (t' 2 ) — (t'} 2 = rt. 

When r — > 0, P(t,t') — > 6(t — t') so that p(t) = p(t) and Vl(<) = exp{— iLi} is the usual unitary evolution. 
However, for finite t, the evolution operator Vr,(i) describes a decay of the off diagonal matrix elements in the energy 
representation, whereas the diagonal matrix elements remain constants, i.e. the energy is still a constant of motion. 
In fact, in the energy eigenbasis, Eqs.d|) and (|]) yield 

P«,m(*) = e-' r ""-*e- i *'-»Vm(0), (6) 

where 

7n,m = — log (l + w^ m r) , (7) 

i/ n m = - arctan (w„. m r) , (8) 

T 

with TiLOn t m — (E n — E m ) the energy difference. One can recognize in Eq.(^), beside the exponential decay, a frequency 
shift of every oscillating term. 

The phase diffusion aspect of the present approach can also be seen in the evolution equation for the averaged 
density matrix p(t). Indeed, by differentiating with respect to time Eq.(||) and using (||) one gets the following master 
equation for p(t) 

-^p(t) = --log(l+iLr)p(t). (9) 
at t 

It is worth noting that by expanding the logarithm up to second order in r, one obtains 



;P{t) = - « [H.TKf)] - ^tt W, W ,?(»)]] . (10) 



which is the well known phase-destroying master equation Hence, Eq.(^) appears as a generalized phase destroying 
master equation taking into account higher order terms in r. Nonetheless, the present approach is different from the 
usual master equation approach, in the sense that it is model independent, non perturbative and without specific 
statistical assumptions. In fact, the probability distribution (|^), is derived only from the semigroup condition, and it 
is interesting to note that this condition yields a gaussian probability distribution as a limiting case when t>r. 

III. ATOMIC DIFFRACTION BY STANDING LIGHT WAVE 

We consider the deflection of an atomic beam of two-level atoms by a classical standing wave in a cavity as in 
Fig. (|l|) . The atomic beam is normally incident on the standing wave and experiences an exchange of momentum with 
the photons in the light wave. We shall assume that the frequency of the light field is well detuned from the atomic 
resonance so that we may neglect the spontaneous emission. The Hamiltonian describing the system is [H 

p2 

H = hwoaz + ^ + m {o-z~' lUt + (T+e luJt ) cos(fcX) , (11) 
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where P = (Px,Py-,Pz) is the center of mass momentum of the atom, M is the atomic mass, az and a± are the 
pseudo spin operators for the atom, loq and u> are the atomic and field frequencies, k — lu/c is the wave number of the 
standing wave, and fl the Rabi frequency. 
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FIG. 1. Scheme 



We shall assume that the interaction time is sufficiently small that the transverse kinetic energy absorbed by the 
atom during the interaction can be neglected (Raman-Nath approximation). This is equivalent to neglect the term 
P 2 /2M in the Hamiltonian (|IT|). Then, transforming to the interaction picture (with Ho ~ Tuooz) the Hamiltonian 
may be written in the form 



H = hAaz + 2htta x cos(kX) , 

= hV(X) [cos 6{X)a z + sin 6{X)a x ] 



(12) 
(13) 



where A = lo - w, V(X) = ^/A 2 +An 2 cos 2 (fcX), cos6(X) = A/V(X) and sm0(X) = 2Clcos(kX)/V(X), 
In the limit of large detuning, i.e. A ^> 2ilcos(kX), the effective Hamiltonian remains 



H = hAa z 
The atomic state vector can be written as 



2hti i cos 2 (kX) 
A 



<JZ 



(X\*(t)) = a(X,t)\e) + b(X,t)\g), 

where \g), and |e) are the ground and excited internal atomic states, while a, and b are probability amplitudes. 
We assume the atoms initially in their ground state with a Gaussian wavefunction 



a(x,0) =0, 

b(x,0) = ( 7 ro- 2 )- 1 / 4 exp 



x 
~2Tx 



(14) 
(15) 

(16) 
(17) 



where x = kX and e x is the transverse position spread of the input beam. 
The Schrodinger equations are 



dt 
dt 



A 
~2 



A 
~2 



Q cos x 
A 

n 2 cos 2 x 



b, 



(18) 
(19) 



Of course, it results a(x, t) = 0, while the solution for b can be written in terms of the Bessel functions of the first 
kind J n (Tol, as 
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b(x, t) = exp 



n 2 \ ' 

2KY 



-Hoc 



*" J « 
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n 2 t\ 



2A J 



e 2mx b{x, 0) . 



(20) 



Taking the Fourier transform of this relationship shows the effect in the momentum space 

+°° / 2^ 



b(p x ,t) = exp 



A W_ 
2~ + 2A 



inj n 



n— — oo 
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2A 



b{ Px -2n,0). 



(21) 



where b denotes the Fourier transform of b, and p x = Px/hk. The scattered ground state wavefunction is a superpo- 
sition of Gaussian modulated plane waves with momentum 2nhk. The momentum transferred from the field to the 
atom is an even multiplies of tik corresponding to the absorption of a photon from the (+k) component followed by 
induced emission in the (— k) component of the standing wave. 

The final output momentum probability, which can be measured, is composed of a comb of images of the initial 
momentum distribution 



-oo +0O 



j(p x ,T)= b(p x ,T) = i {n - m) Jn(T/4)J m (T/<i)b(p x -2n,0)b(p x -2m,0), 



(22) 



n— — oo m— — oo 



where T = 2fl 2 t/A. In order to resolve the peaks it is necessary to have a narrow initial momentum spread, i.e. 
e x > 1. 



IV. EFFECTS OF TIME FLUCTUATIONS 



The atomic wave packet has a finite position spread also in the longitudinal direction (Z in Fig.l). For the sake of 
simplicity we assume e z — e x — e = kez = kex ■ The non-zero value of e 2 give rise to an uncertainty in the arrival 
time H of the atoms at the end of the cavity, i.e. at the time T of Eq.(|2^). We can account for this uncertainty by 
considering the time as a random variable and use the quantum mechanical consistent approach described in Section 
II Q. Then, accordingly to Eq.(|l|), we have to do the following replacement 

w{ Px ,T) ee Tv{ PxP (T)} ~^w(p Xl T) ee Tv{ Px p(T)} . (23) 

From Eqs.© and (|), it results 

+ oo -j-oo 

w( Px ,T)= ]T i {n - m) In, m (T)b(p x -2n,0)b( Px -2m,0), (24) 

n— — oo m— — oo 

where 

r°° i p-t'/t 

In, m {T) = J o dT'j n {T'/±)J m {T'/l)-— r —{T'/T) T ' r -\ (25) 

with T = 2J7 2 r/A, and the integral given explicitly in the appendix. 

The output momentum distribution is shown in Fig. 2 for several values of T. Notice, that the chosen interaction 
time T implies the deflection of all atoms in the ideal situation a), where no peak is present at p x = 0. However, 
the effect of intrinsic decoherence is to reduce the exchange of momentum between atoms and field, then to leave the 
atoms undeflected. In fact the central peak increases from b) to c) while the others decrease. 
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FIG. 2. Momentum distribution w for T = a); T = 1 b); and T = 10 c). The values of other parameters are: T = 10, and 

€ = 10. 



V. CONCLUSIONS 



We know that T is related to the time parameter r. The latter represents the time uncertainty in the interaction 
process. On the other hand, the uncertainty in localizing temporally a wave packet at a given space point, which give 
rise to the time energy uncertainty relation H, may be considered as SZ/v, where 5Z is the width of the packet and 
v its group velocity q|. Then, we may assume H 

where (Pz)/M is the mean velocity of the atoms along the Z direction, while the square root term represents the 
width of the packet along the Z direction. The first two terms inside the square root of Eq.(|2^) describe the well 
known free particle wave packet spread (during the interaction time t), while the last term is due to the classical 
momentum spread. In this case r (hence T) becomes time dependent, but the positive map described by Eqs.(|l|), (||), 
can be maintained by dropping the semigroup property. 

Of course, if in the above expression (^6|) £p z is the dominant term, we only have classical decoherence effect, 
essentially due to the thermal spread of momentum. However, even if such effect is eliminated, decoherence will be 
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still present having a purely quantum nature. Moreover, it will be intrinsic to the system, and not due to external 
degrees of freedom. 

Effects similar to those of Fig. (0) can be recognized in the experimental data of Ref. S , but mainly due to (classical) 
velocity spread. Nevertheless, more recently, cooling techniques offer the possibility to have almost monochromatic 
atomic beam (e.g. 5{P z )/{Pz) = 10~ 3 with {P z )/M = 10 3 ms" 1 @). Then, it would be possible to achieve 
minimum uncertainty states for the atomic wave packet, with e.g. ez ~ 10~ n m Jf2[ . This means that the dominant 
term in Eq.(|26|), once fixed the interaction time t w 1CP 9 s ||, becomes the Schrodinger spread, i.e. the second term 
in the square root. That gives T«l which is promising to test the intrinsic quantum nature of decoherence as can 
be seen from Fig.(||). On the other hand, the opposite limit, where ez is dominant in Eq.(|2^) could be investigated 
as well with slow atoms at large de Broglie wavelength [|l3| . 

In conclusion, we have studied intrinsic decoherence arising from the interaction time fluctuations in atomic diffrac- 
tion by standing light wave. The main limitation to observe such effect relies on the necessity to suppress classical 
fluctuations, however, the technological achievements in atom optics, atom lasers |l4| at last, make the possibility to 
observe intrinsic quantum decoherence realistic. Finally, it is worth interesting to also apply the above theory to the 
problem of atomic wave diffraction in time fl5[] ; this is planned for future work. 



APPENDIX 

The integral (^5|) can be written explicitly as follows 

In,m(T) = I\n\,\m\(T) (_)(l»|-n+l-l— )/2 , ( 27 ) 

where 

T(M+M)r (|„| + |m| + $ ) 4 F 3 (u, v, -£) 

2 3(M+M) r ( 1 + |„|)r(i + | m |)r (f) 

with F the generalized hypergeometric function Jl(]] whose vectors u and v are 



I\n\,\m\(T) - o3flnl + lmnTVl , I^IYrVl , \^\\T(T\ ' ( 28 ) 



_ (I \n\ \m\ 1 \n\ \m\ \n\ \m\ T 1 \n\ \m\ T \ 
v = (1 + \n\, 1 + H, 1 + |n| + \m\) . (30) 
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